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The emission from an electron in the field of a relativistically strong laser pulse is analyzed. At 
pulse intensities of J > 2 • 10 22 W/cm 2 the emission from counter-propagating electrons is modified 
by the effects of Quantum ElectroDynamics (QED), as long as the electron energy is sufficiently 
high: £ > 1 GeV. The radiation force experienced by an electron is for the first time derived from 
the QED principles and its applicability range is extended towards the QED-strong fields. 

PACS numbers: 52.38.-r Laser-plasma interactions, 41.60.-m Radiation by moving charges, 52.38.Ph X-ray, 
gamma-ray, and particle generation 
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I. INTRODUCTION 



A. QED strong fields 



In Quantum ElectroDynamics (QED) an electric field, 
E, should be treated as strong if it exceeds the Schwinger 
limit: E > E$, where 



E S = 



m e c 



HA 



c 



(see [l|). Such field is potentially capable of separating a 
virtual electron-positron pair providing an energy, which 
exceeds the electron rest mass energy, m e c 2 , to a charge, 
e = — |e|, over an acceleration length as small as the 
Compton wavelength, 

A c = — w 3.9 ■ Kr n cm. 

m e c 

Spatial scales associated with the field should be greater 
than Ac- 

Typical effects in QED strong fields are: electron- 
positron pair creation from high-energy photons, high- 
energy photon emission from electrons or positrons and 
electron-positron cascade development (see Q- Q) as the 
result of the first two processes. 

Less typical and often forbidden by conservation laws 
is direct pair separation from vacuum. This effect may 
only be significant if the field invariants as defined in [J] , 
Fi = (B • E), F 2 — E 2 — B 2 , are large enough. Indeed, 
the considerations relating to pair creation are applicable 
only in the frame of reference in which B = or B||E. 
The electric field in this frame of reference, Eq = F 2 /2 + 
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a/F 2 + F 2 /4, exceeds the Schwinger limit only if the field 
invariants are sufficiently large. 

Here the case of weak field invariants is considered: 



\F, 



< E 2 



(1) 



and any corrections of the order of Fi/Eg, FijE\ are 
neglected (see [5] about such corrections). So, neither 
the cases when the field itself is too strong, nor the cases 
when its spatial scale is too short are considered here. 
Below, the term 'strong field' is only applied to the field 
experienced by a particle (electron or positron). 
Particularly, a QED-strong electric field, 



En 



|pxB| 

m e c 



(2) 



may be exerted on relativistic charged particles with mo- 
mentum, p, gyrating in the strong magnetic field, B, of a 
neutron star, as the result of the Lorentz transformation 
of the electromagnetic field. The field as in Eq.® may 
exceed the Schwinger limit, as long as |p| 3> m e c and/or 
the magnetic field is strong enough. 



B. QED-strong laser fields 

In a laboratory experiment QED-strong fields may be 
created in the focus of an ultra-bright laser. Consider 
QED-effects in a relativistically strong pulsed field [||: 



/ a 2 > 1, 



(3) 



A being the vector potential of the wave. In the labo- 
ratory frame of reference the electric field is not QED- 
strong for achieved laser intensities, J ~ 10 22 W/cm 2 [y], 
and even for the J ~ 10 25 W/cm 2 intensity projected 
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Q. Moreover, both field invariants vanish for ID waves, 
reducing the probability of direct pair creation from vac- 
uum by virtue of the laser field's proximity to ID wave. 

Nonetheless, a counter-propagating particle in a ID 
wave, a(£), £ = u)t — (k • x), may experience a QED- 
strong field, Eq = |dA/<i£;|w(£ — p\\)/c, because the laser 
frequency, uj, is Doppler upshifted in the frame of refer- 
ence co-moving with the electron. Herewith the electron 
dimensionless energy, £ , and its momentum are related to 
m e c 2 , and m e c correspondingly, and subscript || herewith 
denotes the vector projection on the direction of the wave 
propagation. The Lorentz-transformed field exceeds the 
Schwinger limit, if 



X ~ Eq/Es 



Ac 
A 



(£-P||) 



da 



> i, 



(4) 



where A — c/uj. Note, that the above-mentioned restric- 
tion on the field spatial scale is here assumed to be ful- 
filled for the upshifted wave frequency: 



cj(£-pu) < c/X c . 



(5) 



Nevertheless, the condition as in (0| may be fulfilled as 
long as the field is strong enough. Numerical values of the 
parameter, x, may be conveniently expressed in terms of 
the local instantaneous (not time-average!) intensity of 
the laser wave, J: 



X 



- — (£ 

2 A [t 



P\\ 



da 



0.7 



J 



(£-pw) r 

10 3 V 10 23 [W/c 



the choice of the numerical factor of 3/2 is explained be- 
low. For a counter-propagating electron of energy ~ 1 
GeV, that is for (£ -py) ~ 4- 10 3 , the QED-strength pa- 
rameter is greater than one even with the laser intensities 
already achieved. 

The condition of x > 1 also separates the parameter 
range of the Compton effect from that of the Thomson 
effect, under the condition of Eq.([3]). The distinctive fea- 
ture of the Compton effect is an electron recoil, which 
is significant, if a typical emitted photon energy, Tluj c , 
is comparable with the electron energy Q. Their ratio, 
X = Xclo c / (c£), equals x as defined in Eq. ([7]) with the 
proper numerical factors (cf Eq. lfTS]) ). It should be noted, 
however, that under the conditions discussed in (|3l5j) the 
Compton effect drastically changes. 



C. Classical radiation loss rate as an input 
parameter for QED 

Radiation processes in QED-strong fields are entirely 
controlled by the local value of Eq (this statement may 
be found in §101). A good signature for Eq is the ra- 
diation loss rate of a charge, as introduced in the classical 
electrodynamics: 



I c l(E ) 



2e A El 
3m 2 c 3 



3m 2 c 3 



(6) 



which is a Lorentz invariant. Therefore, it may be ex- 
pressed in any frame of reference in terms of a 4-square 
of the Lorentz 4-force, p = (/°,f) = £{i^ ■ v/c,fW), 
where v is the velocity vector and: 



f (3) = eE + - [v x B] 



/° = eE-p, f = e£E + e[p x B]). 

So, the QED-strength of the field may be determined in 
evaluating I c \ and its ratio to Ic = 1 C i(2-E g/3): 



X 



Ic 



ie 2 c 



27Xi 



(7) 



If X > 1 then the actual radiation loss rate differs from 
7 c l, however, it may be re-calculated using I c \ as a sole 
input parameter. 



D. Possible realization of QED-strong fields in 
laboratory experiments 

The first experiments which demonstrated QED ef- 
fects in a laser field were fulfilled using an electron beam 
of energy « 46.6 GeV (see [9(), which interacted with 
a counter-propagating terawatt laser pulse of intensity 
J ~ 10 18 W/cm . A reasonably high value of x ~ 0.4 
had been achieved, however, the laser field was not rel- 
ativistically strong with |a| < 1. The high value of x 
had been achieved at the cost of very high energy of 
the upshifted laser wave: the transformed photon energy 
amounted to ~ 0.1 MeV, which is not small as compared 
to the electron rest mass energy, m e c 2 f=s 0.51 MeV. It 
could be interesting to upgrade this experiment towards 
the highest achievable laser intensities > 2 • 10 22 with 
the use of a wakefield-accelerated beam of electrons of 
energy ~ 1 GeV (see [HI). First, the 2-3 times larger 
value of the QED parameter, x-i may be achieved with 
the exponentially increased probability for pair creation. 
Second, such experiment would be highly relevant to the 
processes which will occur in the course of laser-plasma 
interaction at even stronger laser intensities. 

Indeed, counter-propagating electrons can be gener- 
ated while a laser pulse is interacting with a solid tar- 
get. For this reason, the radiation effects in the course 
of laser-plasma interaction are widely investigated ((see 
8, 12]). With future progress in laser technology and 
by achieving intensities of J ~ 10 23 W/cm laser-plasma 
interactions will be strongly modified by QED effects, so 
that the capability to model these effects now is of inter- 
est. 



E. Radiation back-reaction 

The principle matter in this paper is an account of the 
radiation back-reaction acting on a charged particle. The 
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radiation losses reduce the particle energy, affecting both 
the particle motion and the radiation losses themselves. 

This effect can be consistently described by solving 
the dynamical equation which appears to reduce to the 
modified Lorentz-Abraham-Dirac equation as derived in 
13]. 14]. The new element here is that in this equation 
the radiation back-reaction on the electron motion should 
be expressed in terms of the emission probability, while 
applied to QED-strong fields. 

In section II the emission from an electron in the field 
of a relativistically strong wave is discussed within the 
framework of classical electrodynamics. The transition 
from the vector amplitude of emission as described in 
[TEj to the instantaneous spectrum of emission is treated 
in terms of the formation time, a concept, which is not 
often used in classical electrodynamics. For QED-strong 
laser pulses the calculation of the emission probability 
is given in Section III. The radiation processes in QED- 
strong fields appear to be reducible to a frequency down- 
shift in the classical vector amplitude of emission result- 
ing from the electron recoil, accompanied by a contribu- 
tion to emission associated with the magnetic moment of 
electron. The radiation effect on the electron motion in 
strong fields is discussed in Section IV. The Conclusion 
summarizes the results and discusses future prospectives. 



A. Transformed space-time 

A method facilitating many derivations involves the 
introduction of a specific time-space coordinate frame. 
Consider a ID wave field taken in the Lorentz calibration: 

a" = a M (0, 6 = (* ' x), (k-a) = 0, 

a M = (0,a), fc M and x^ being the 4- vectors of the poten- 
tial, the wave and the coordinates. Herewith the 4-dot- 
product is introduced in a usual manner: 

(k ■ x) — k^x^ = urt — (k • x) 

etc. Space-like 3- vectors (i.e., the first to the third com- 
ponents of a 4-vector) in contrast with 4-vectors are de- 
noted in bold, 4-indices are denoted with Greek letters. 
Note that a metric signature (+,—,—,—) is used, there- 
fore, for space-like vectors the 3D scalar product and 4- 
dot-product have opposite signs, particularly: 



/da 



Introduce a Transformed Space-Time (TST) : 



(cf =F x\\ 



2, 



„2,3 



X_L, 



II. EMISSION IN RELATIVISTICALLY 
STRONG FIELDS 



In this section QED effects are not yet considered, but 
the electromagnetic wave field is assumed to be relativis- 
tically strong. Angular and frequency distributions of 
electron emission are discussed. The goal is to estab- 
lish a connection between the methods usually applied 
to calculate emission in weaker fields, on one hand, and 
the conceptually different QED approach on the other. 
For relativistically strong laser fields, even though QED 
effects do not yet come into a power, still some concepts 
of the QED emission theory appear to be applicable and 
useful, among them are the formation time of emission 
and instantaneous spectrum of emission. 

In weaker fields, especially for the particular case of 
a harmonic wave, the emitted power is given by an in- 
tegral over many periods of the wave. This standard 
approach, however, may become meaningless as applied 
to the ultra-strong laser pulses, for many reasons. These 
pulses may be so short that they cannot be thought of 
as harmonic waves. Their fields may be strong enough 
to force an electron to expend its energy on radiation 
faster than a single wave period. However, an even more 
important point is that the radiation loss rate and even 
the spectrum of radiation is no longer an integral charac- 
teristic of the particle motion through a number of wave 
periods: a local dependence of emission on both particle 
and field characteristics is typical for the strong fields. 



subscript _L denoting the vector components orthogonal 
to k. The properties of the TST provide a convenient 
description for the classical motion of an electron in the 
ID wave field. First, note, that 



V2' 



P 



Hk-P) 
V2 ■ 



(p-k) 



Second, the generalized momentum components, p° and 
pxo = P_l + a, are conserved. Third, the metric tensor 
in the TST is: 

G01 s~<10 i y^f22 ^33 i /^M^ 

— ' i — 1, u — u — 1, v ' — ^ fiv • 

1, being expanded in the 



Finally, the identity, £ 2 = p 2 
TST metric, gives: 



1 



P 



pi 



2p° 



1 + (P-LO ~ a) 
V2X{k -p) 



The classical radiation loss rate is found by virtue of ex- 
panding the Lorentz force squared in the TST: 



/rJ = 



2e 2 (/■/) = 2e 2 c{k-pf 
3c m 2 c 2 3 



da 

di 



(8) 



The derivative over x° or, the same, over £ is conveniently 
related to the derivative over the proper time for electron: 



dr 



= £ 



a . d ; 
dt + {v -d^ ] 



= c(k ■ p) 



(9) 
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B. Classical trajectory and momenta retarded 
product 

Many characteristics of emission may be expressed in 
terms of the relationship between the 4-momenta of the 
electron at different instants: 

p"(o =p^') - sat* + 2(p(r) ' ;f a) : (fo) V, (io) 



2{k-p) 



5a» = a^(0 - a»(?). 



where 



As a consequence from Eq. (fl0|) . one can obtain the ex- 
pression for the Momenta Retarded Product (MRP): 



(p(0-pO = 1-^ = 1 



(6a) 



(11) 



Note, that the MRP is given by Eq. flTTj) for an arbi- 
trary difference between £ and but only for the par- 
ticular case of the ID wave field. However the limit of 
this formula as |£ — £'| — > 0, which is as follows: 



1 



'\2 



(KO-p(0)lis-mo«i + 2«-£') 



da 



or, in terms of the MRP in the proper time, r: 



(p(r) -p(t + St)) = 1- ($r) 



a (/•/) _ 
2m?c 2 



(12) 



has a much wider range of applicability. Eq. (fT2f is de- 
rived from the equation of motion: 

dp» _ /" 
dr m e c 

using the identities: 

( P (t)-p(t)) = 1, (p(r) • /(r)) = 0, 



d(p(r) ■ p(r + <5r)) / dp dp\ 2 

— — = ^ T [d? ■ ^ J + ° {{5t) ] - 



Here the superscript asterisk means the complex conju- 
gation and the vector amplitude of emission, A c i(u/), is 
given by the following equation: 



A c i(u/, n) = - / v(i) exp < iui'[t 

C J — I C 



(n-r(t)). 



df. 



see Eq. (14.67) in Ref.[15[ followed by the discussion of the 
way to account for a polarization. The use of the same 
notation, A, both for the emission vector amplitude and 
for the vector potential should not mislead the reader. 
Recall, that the emission vector amplitude is closely re- 
lated to the Fourier-transformed vector potential in the 
far-field zone of emission. 

Introduce a 4- vector amplitude of emission, 



A>',n)=e 



p^ir) exp 



ic / (A/ • p{r'))d7 



dr, 



which is expressed in terms of the proper time for the 
electron and its (dimensionless) 4-momentum. As long 
as (fc' • p)/c is a frequency of the emitted photon in the 
frame of reference co-moving with the electron, the 4- 
vector amplitude is the Fourier integral of the electron 
4-momentum with the Lorentz-modified frequency. 

Note the following properties of the 4-vector ampli- 
tude. First, its space-like vector components, which are 
perpendicular to the wave vector of the emitted pho- 
ton, coincide with those for 3-vector amplitude, hence, 
they quantify the polarization properties of the emission 
for two different polarizations. Second, the dot-product, 
(A c \ ■ k'), vanishes as being the integral of a perfect time 
derivative. Now construct the dot-product, (A c \ ■ A*J 
and expand it in the TST, which is formulated in terms 
of the emitted wave: 



(A c i ■ A*j) — A cl (A cl )* + Al^A^)* 



\A 2 cl \ 2 



\A 3 cl \ 2 . 



From the above properties of the 4-vector amplitude, the 
first two terms vanish identically as A^ oc (A c \ ■ k') = 
0. The other two terms as taken with the proper factor 
give the emitted energy summed up over polarizations, 
therefore, the latter sum may be expressed as follows: 



E 



dR 



cl 



doj'dn 



4ir 2 c 



(Ad ■ A* cl ). 



(14) 



C. Vector amplitude of emission in classical 
electrodynamics 

In Ref.[l5| the frequency spectrum and angular dis- 
tribution, dR c \/ (duj' dn), of the radiation energy, dR c \, 
emitted by an electron and related to the interval of fre- 
quency, du)', and to the element of solid angle, dn, for 
a polarization vector, 1, is described with the following 
formula: 



di? c 



duj'dn Air 2 c 



U ' } |(A cl (o;')-r)| 2 . 



(13) 



Now introduce the radiation loss rate, dl c i/(du>'dn) re- 
lated to the unit of time, the element of a solid angle and 
the frequency interval. Its connection to dR c \ / '(duj 1 'dn) is 
evident: 



E 



dR, 



cl 



dw'dn 



+oo 



dldjt) 
dw'dn 



dt = 



+ 00 



dI c x{T) 
dui'dn 



£(r)dT. 



In Eq.([T4| the dot- product of the 4-vector amplitudes, 
(A c \ ■ A*j) is in fact the product of two integrals over dr, 
which can be represented as the double integral, over, 
say, dr\dT2- Transform the integration variables in this 
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double integral, by introducing r = (n +r 2 )/2, = n — 
t 2 . The spectral and angular distribution of the radiation 
loss rate may be expressed in terms of the Fourier integral 
of the MRP: 



jjdW = e 2 (a/) 2 
rfw'dn 47r 2 cf(r) 



+ -) -p(r- -))x 



x exp 



r+e/2 

ic I (k' ■ p(r'))dr' 

r- 0/2 



(20. 



D. Frequency spectrum and formation time 

The specific feature of the particle relativistic motion 
in strong laser fields is that the main contribution to the 
above integral comes from a brief time interval with small 
values of 0. A closely related point is that the emitted 
radiation is abruptly beamed about the direction of the 
velocity vector, p(r)/|p(r)|. Therefore, in the following 
expansion of the frequency in the co-moving frame, 



(k'.p(r')) = [k' 



-^p(r)]-p(r')j+^(p(r)-p(r')) 



in the first dot-product one may approximate p M (r') w 
P^(t). Indeed, for an ultrarelativistic electron and a pho- 
ton, which both propagate in the same direction, the 
difference between p^/£ and c(fc') M /u/ is already small, 
therefore in the second multiplier of the dot-product the 
small difference, p M (r) — p^(r'), may be neglected: 



h0/2 

c / (k' ■ p(t'))<1t' « 9c(k' - P (t))- 

It-9/2 



he/2 



-0/2 



[(p(r)-p(r'))-l]dr'. 



Now the only angle-dependent multiplier is exp[— i6c(k' ■ 
p{r))]. For simplicity, the angular spectrum of emission 
can be approximated with the Dirac function: 

(2/ c i(r) _ $2 f n _ J>\ dldir) 



duj'dn \ |p| / dui' 
and with the use of the formula (see §90 in Ref . ( fiol] ) . 



x sm ■ 



S(t) 



r+e/2 



T-0/2 



[(p(T)-p(r'))-l}dr' 



(20. 



Thus, the frequency spectrum of emission is entirely de- 
termined by the MRP, which is a scalar Lorentz-invariant 
funcion of the proper time. Both the fore-exponential 
factor and the argument of the exponential function de- 
pend on the mentioned MRP. Therefore, both the spec- 
tral composition of the MRP and its magnitude may be 
of importance. Their relative role is controlled by the 
ratio of the frequency of the electron motion, ojq, to the 
acceleration magnitude, both being determined in the co- 
moving frame of reference. Here the field is assumed to 
be so strong, that the acceleration it causes plays the 
dominant role, i.e. the following inequality is claimed: 



(/•/) 



» < 



(15) 



Under these circumstances, the integral determining the 
emission spectrum is calculated by virtue of the displace- 
ment of the integration contour in the plane of the com- 
plex variable, 0, so that the deformed contour passes 
through the point of a stationary phase, s t ■ In this 'sad- 
dle' point the phase gradient turns to zero: 



rf0 



•+0/2 



-0/2 



[(p(T)-p(r'))-l}dr' 



The larger the acceleration becomes, the closer the sta- 
tionary phase point, st draws to the real axis, and, 
hence, the shorter the time interval becomes, ~ 0f = 
|0 s t|, which gives the non- vanishing contribution to the 
emission spectrum. The characteristic duration of this 
time interval, 0f = |0 s t| is referred to as a formation time 
(or coherence time - see [H|],[l3])- At the limit of large 
accelerations the formation time is given by the following 
formula: 



±i 



2m e c 



2m e c 



vM/ 7 /) 



(16) 



where the approximation for the MRP as in Eq. (fT2j) is 
applied at |0| < 0f . With the use of Eq.(fT12j) the universal 
emission spectrum is obtained: 



(2/ cl (r) e 2 to' 
du? ~ 2ttc£ 2 (t) 



+00 



1 



(/(t)-/(t))0 1 



2m 



2„2 



exp[i0c(fc' • p(r))]dn — y^x7> ex P ' 



2S(t) 



the following expression may be found for the frequency 
spectrum of emission: 



rfJci(r) 



2ttc£ 2 (t) 



+ 00 



\(p{t + 5) -P( T - 2)) x 



£{r) 



(/(r)-/(r))0 3 
24m 2 c 2 



(20. 



The integral can be expressed in terms of the MacDonald 
function (= the modified Bcssel function): 



dUj) = Ici(t)_ 
duo' uj c 



Q c \(ro), 



dro 



Id{r)Qd{ro), 
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where Q c i{fo) is the unity- normalized spectrum of the 
gyrosynchrotron emission (J Q c \(r)dr = 1): 

th/3 r°° u)' 
Qcl(ro) = ^r / K 5/3 (r')dr', r = -, (17) 

and 

3 g(T)V-(/(r)-/(r)) 3 ,,/ 3/ cl (r)c 
" c = 2 ^ = 2 £(r) V —2^—' 

or, which is the same, 

^=eM±=e x . (is) 

Note, that despite all approximations, the integral over 
the frequency spectrum is consistently equal to I c \. 



is fine structure constant. The ratio (12~T1) is much less 
than unity as long as % < 1. Note that in the oppo- 
site limiting case of QED-strong field, the extra factor 
of Jqed/^cI oc x ( see FigEl below) makes ratio (|2"Tj) 
small at x ^ 1 as well. 

The same estimate of the formation time is applicable 
to any relativistically strong electromagnetic field, not 
only to ID wave. Particularly, in the wakefield accel- 
eration scheme, where the electric field in the wakefield 
of the pulse may be even larger than the relativistically 
strong field in the pulse itself. With this account, the ac- 
celeration of almost monoenergetic electron beam by the 
laser pulse must be accompanied by the gyrosynchrotron- 
like spectrum of emission (which are actually observed - 
see [l8|, [Hj]). These observations demonstrate the gen- 
eral character of the gyrosynchrotron emission spectrum 
(this point of view, presented in [18| , may be also found 
in §77 in [4]). 



E. Implications for strong laser fields 

As discussed, the condition £ > 1 and the inequality 
(fl~5|) are both fulfilled for an ultra-relativistic electron 
gyrating in a uniform steady-state magnetic field. By 
expressing the 4-force squared, — (/ • /) = e 2 p^B 2 (see 
Eq-©), and taking the gyrofrequency in the co- moving 
frame, w 2 = £ 2 e 2 B 2 /(?n 2 c 2 £ 2 ) = e 2 B 2 /(m 2 c 2 ), one finds 
that (fT5|) is fulfilled as long as 3> 1. 

Furthermore, application to the ID wave field is no less 
straight forward. The laser wave frequency in the comov- 
ing frame, ujq = c(k ■ p), is present on the left hand side 
of Eq.©. The Lorentz 4-force squared is given in Eq.©, 
resulting in the following estimate for the formation time: 



9{ = 



c(k-p)\d a /dt\ 



(19) 



Now it is easy to see, that the condition, 8[UJo <C 1 as in 
Eq. (fl~5|) is fulfilled in relativistically strong wave field, at 



da 



> i. 



(20) 



The formation time tends to zero as the wave ampli- 
tude tends to infinity. The change in the electron energy 
within the formation time is always much less than the 
particle energy, £m e c 2 . Within the classical field theory 
this statement follows from Eqs. (|8ll9[) . With an account 
of an extra factor of £ , which arises while transforming 
the formation time to the laboratory frame of reference, 
the relative change in energy equals: 



m e c 2 



where 



ot 



e 

Tic 



1 
137 



(21) 



Emission within short time intervals and 
implications for numerical simulation 



In strong fields satisfying the condition as in (|T5|) both 
emission vector amplitude and the emission spectrum 
may be determined with respect to brief time interval, 
At, which may be much shorter than the field period. 
The only requirement is that this interval should be large 
as compared to the formation time: 



t+At 



df 

W: 



»0 f , 



(22) 



however, the change in the field and particle character- 
istics within this time interval may be small, as long as 
Atu <C 1. In this case the span in the integral determin- 
ing the vector amplitude may be chosen to be (t, t + At). 
However, in the integral over d9, which determines the 
emission spectrum, the integration limits are much larger 
than the formation time, therefore, they may be again set 
to (— oo, +00). 

These considerations justify the numerical scheme for 
collecting the high frequency emission as described in 
[lij (which does not seem different from that briefly de- 
scribed in 18]). In addition to calculating the electro- 
magnetic fields on the grid using Particle-In-Cell (PIC) 
scheme, in which fields are created by a moving particle 
within the time step, At, one can also account for the 
higher-frequency (subgrid) emission spectrum, by calcu- 
lating the instantly radiated energy, I c \At, and its dis- 
tribution over frequency, parameterized via I c \. Another 
often used approach based on the calculation of the vec- 
tor amplitude of emission (see, e.g.,[2(|) seems to be less 
efficient, although, theoretically, should provide the same 
result. The vector amplitude formalism, on the other 
hand, may be better applicable to the cases, where the 
high-fre quen cy emission from multiple electrons is coher- 
ent (see 21|). Stemming from these considerations it is 
now easy to proceed to the QED approach. 
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III. ELECTRON IN QED-STRONG FIELD: THE 
EMISSION PROBABILITY 



The emission probability in the QED-strong ID wave 
field may be found in §§40,90,101 in [lGj, as well as in 
16], 22]. In application to the wakeficld acceleration of 
electrons of energy 1 TeV the QED effects had been 
also discussed in [23]. However, to simulate highly dy- 
namical effects in pulsed fields, one needs a reformulated 
emission probability, related to short time intervals (not 
(—oo, +oo)). 

Indeed, it is demonstrated above that in strong fields 
the emission processes are essentially local functions of 
the instantaneous parameters. Therefore, in QED-strong 
fields the emission probability should be formulated in 
terms of the local values of the electromagnetic field in- 
tensities, or, the way we adopt, it may be parameterized 
via the classical radiation loss rate or the Lorentz 4-forcc 
squared: — oc I c \. This emission probability is red- 
erived here with careful attention to consistent problem 
formulation and neglecting technical details. 



A. A QED solution of the Dirac equation 

The Dirac equation which determines the evolution of 
the wave function, for a non- emitting electron in the 
external field, reads: 



(23) 



7' 



being the Dirac 4*4 matrices, (7 , 7 1 , j 2 , j 3 ). The rel- 
ativistic dot-product of the Dirac matrices by 4-vector, 
such as (7-0), is the linear combination of the Dirac 
matrices: (7-0) = j°a° — 7 1 a 1 — j 2 a 2 — j 3 a 3 . Such 
linear combination, which is also a 4 * 4 matrix, may 
be multiplied by another matrix of this kind or by 4- 
componcnt bi-spinor, such as ip, following matrix multi- 
plication rules. For example, (7 • ajtp is a bi-spinor, as is 
the matrix, (7 • a) multiplied from the right hand side by 
the bi-spinor, ip. 

By expanding Eq.Q23j) in the TST it is easy to find its 
solution in a form of a plane electron wave (the normal- 
ization coefficient N — const): 



VN 



V2 



(24) 



Here u(p(£)) is plane wave bi-spinor amplitude, which 
satisfies the system of four linear algebraic equations: 



(7-p(0Mp(0) = «(p(0)> 



(25) 



as well as the normalization condition: iiu — 2. The .in- 
dependent momentum, f>(£), in the bi-spinor amplitude 
should be taken in accordance with Eq. fTU]) as for the 



classical trajectory of the electron, 
phase multiplier, P(£) is as follows: 



The £- dependent 



P(0 = exp - 



P(£) = P(£')exp 



Using Eq. fTO]) . one can find: 



i + pj(£ 2 ; 

2{k-p) 



d£,2 



c i + pj(6; 

, 2{k-p) 



#3 



(26) 



u(p(0) = 



(7-fc) ( 7 - [a(0-oO 



2{k-p) 



u(p(0) (27) 



and verify that Eq. (|2~4"|) satisfies the Dirac equation. The 
advantage of the approach used here as compared to the 
known Volkov solution presented in §40 in [l(| is that 
the wave function in Eqs. (|24H2"T]) is described in a self- 
contained manner within some finite time interval, (£', £) 
(in fact, this interval is assumed to be very short below) 
in terms of the local parameters of the classical trajec- 
tory of electrons. This approach is better applicable to 
strong fields, in which the time interval between subse- 
quent emission occurrences, which destroys the unper- 
turbed wave function, becomes very short. 



B. The matrix element for emission 

The emission problem is formulated in the following 
way. The electron motion in the strong field may be 
thought of as the sequence of short intervals. Within 
each of these intervals the electron follows a piece of a 
classical trajectory, as in Eq. (flO|) . and its wave function 
(an electron state) is given by Eq. (|2"4"l) . The transition 
from one piece of the classical trajectory to another, or, 
the same, from one electron state to another occurs in a 
probabilistic manner. The probability of this transition, 
which is accompanied by a photon emission is calculated 
below using the QED perturbation theory. 

The only difficulty specific to strong pulsed fields is 
that the short piece of the electron trajectory is strictly 
bounded in space and in time, while the QED invari- 
ant perturbation theory is based on the 'matrix element', 
which is the integral over infinite 4-volume. 

To avoid this difficulty the following method is sug- 
gested, which is analogous to the dipole emission theory 
as applied in TST. Introduce domain, A 4 x = (Ax 1 *S±)* 
Ax° , bounded by two hypersurfaces, £ = £_ and £ = 
(see Fig.l). The difference £+ — £_ is bounded as de- 
scribed below, so that A 4 a; covers only a minor part of 
the pulse. A volume, 



d^2 



is a section of A 4 x subtended by a line t — const. 
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FIG. 1: The volume over which to integrate the matrix ele- 
ment while finding the emission probability: in the standard 
scheme for the dipole emission (in dashed lines) and in the 
TST (in solid lines). Arrows show the direction, along which 
the Heisenberg operator advances the wave functions. 



With the following choice for the normalization coeffi- 
cient in Eq.([M|): 



N = 2S±X 



€+ 



the integral of the electron density in the volume V, 

/■«+ . 

i/rfi/idV = S X X ^Vdfc, 



is set to unity, i.e. there is a single electron in the volume 
V. This statement follows from Eq. ([24l) and the known 
property of normalized bi-spinor amplitudes: u ■ 7 • u — 
2£. Here the hat means the Dirac conjugation. 

For a photon of wave vector, (fc') Al , and polarization 
vector, l^, introduce the wave function: 



, A ,y = exp[-i(k' -x)/X c ] p 



or, by expanding this in the TST: 



< A 'Y = 



exp 



where: 



P p (£) = exp 



X C 



iX(k ■ k^x 1 
V2X C 



J", 



2(k-k')X c 



Here the photon momentum and photon energy are re- 
lated to m e c and m e c 2 correspondingly, or, equivalently, 
dimensionless (fc') Al equals dimensional multiplied 
by Xc- The choice of the normalization coefficient, 



N.„ 



2irhcXc 

corresponds to a single photon in the volume, V. 



The emission probability, dW , is given by an integral 
over A x: 



dW = 
Here 



ah fLp 
he 



Li, 



i>f(l ■ (A')*)ip i dx°dx 1 dx 2 dx i 



Vd 3 k' _ hcN p d 2 k' ± d(k ■ k') 



(27rA c ) 3 (27rA c ) 2 (fc • k') 



(28) 



(29) 



is the number of states for the emitted photon. The 
transformation of the phase volume as in Eq. ([29)) is based 
on the following Jacobian: 

/ dk' \ 



d(k' ■ k) 



k' —const 



{k' ■ k)'' 



which is also used below in many places. A subscript i, f 
denotes the electron in the initial (i) or final (f) state. 
The number of electron states in the presence of the wave 
field, Li j, should be integrated over the volume V 



(2ttX c ) 3 



d 3 PlJ dV 



d(k ■ p)ij-d 2 p_ LiJ N i j 
2{2-KfX%{k-p) iJ 



Conservation laws 



The integration by dx 1 dx 2 dx 3 — c\/2dtd 2 x± results in 
three 5— functions, expressing the conservation of totals 
of pj_ and (k ■ p), for particles in initial and final states: 



P± f + k' ± , (k- Pi ) = {k-pf) + {k-k'). 



Twice integrated with respect to da; 1 , the probability dW 
is proportional to a long time interval, At = Ax 1 / (c^/2), 
if the boundary condition for the electron wave at £ = £_ 
is maintained within that long time. On transforming the 
integral over dx° to that over <i£, one can find: 



..d 4 x 



(2TrXc) 3 S ± cAtX 



...de 



x6 2 ( P±l -p± f - k' ± )S((k ■ Pi ) - (k ■ Pf ) - (k ■ k')). 

To take the large value of At seems to be the only way 
to calculate the integral, however, the emission probabil- 
ity calculated in this way relates to multiple electrons in 
the initial state, each of them locating between the wave 
fronts £ = £_ and £ = during much shorter time, 



= (1/c) / £(Od£a/(fr-Pi). (30) 



For a single electron the emission probability becomes: 
dWfi (£_,£+) =StdW/At. 
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Using 5— functions it is easy to integrate Eq. ((28|) over 
dp ±f d(k -p f ): 



dW f 



d(k ■ k')d 2 k' ± (47rA c ) 2 (fc ■ k')(k ■ Pi )(k - Pf ) 
where 



T(0 = 



exp 



Xc(k-pf) 



Pi(£) and P/(£) are the electron phase multipliers, P(£), 
for the electron in initial and final states and 

Prior to discussing Eq.([3"Tj). return to Eq. (j2"T)l and ana- 
lyze it component-by-component in the TST. It appears 
that three of the four components of that equation de- 
scribe the conservation of (k ■ p) and pj_o = p± + a 
for electron in the course of its emission- free motion. 
At the same time, yet another component of Eq. (|27| . 
specifically, p 1 , directed along describes the energy- 
momentum exchange between the electron and the ID 
wave field, maintaining the identity, (p ■ p) — 1. Now 
turn to Eq. (|3"Tj) . Again, three of the four components ex- 
press the conservation of the same variables in the course 
of the photon emission, while the p 1 component, directed 
along fc M describes the absorption of energy and momen- 
tum from the wave field in the course of the photon emis- 
sion. Note, that in the case of a strong field, the energy 
absorbed from field is not an integer number of quanta, 
and that for short non-harmonic field it is not even a 
constant, but a function of the local field. 



D. Calculation of the matrix element 

To calculate the matrix element, one can re- write it 
as the double integral over d^d^i and then reduce the 
matrices u(pij(^)) <S> «(Pi,/(6)) m the integrand to the 
polarization matrices of the electron at £ or at £i using 
Eq. (l2"Tl) . These standard manipulations with the Dirac 
matrices are omitted here. Although in a strong wave 
electrons may be polarized (see [HI), in the present work 
the emission probability is assumed to be averaged over 
the electron initial polarizations and summed over its fi- 
nal polarizations. The ultimate result of these derivations 
is as follows: 



dW 



fi 



, (32) 



d(k • fc')d 2 k^ (27rA c ) 2 0fc • k')(k ■ Pi )(k ■ P} ) 
where 

[(?(£)• p(£i))-i](i-Q<) 2 



D = {l*-Pi{£i)){l-Pi{0) 



AC f 



and 



C 



fi 



(k-Pf) 
(k-Pi) 



1 



Xc(k-k') 
(k ■ Pi) 



< 1 



is a recoil parameter which characterizes the reduction in 
the photon momentum due to emission. 

The matrix element may also be summed, if desired, 
for two possible directions of the polarization vector. The 
second term in the integrand is simply multiplied by two, 
while in the first one the negative of the metric tensor 
should be substituted for the product of the polarization 
vectors (see §8 in [ltj), so that — (pj(£) • f>i(£i)) substi- 
tutes for (I* ■ Pi(£i))(l ■ Pi(£,))- The latter may be trans- 
formed using Eq. dlll) . thus, giving: 



H D = - 



Kp(0-p(£i))-i](i 



2C fl 



E. Vector amplitude of emission in QED case 

Now moving to connection between the obtained re- 
sult, on one hand, and the way the high-frequency emis- 
sion is treated in the framework of classical theory, on 
the other. To facilitate the comparison, both here and in 
Section II the photon frequency and wave vector, u/ and 
k', are not dimensionless. It appears that the QED result 
obtained above can be reformulated in a form similar to 
Eq. (fT3")) . Using the following relationships between the 
differentials: 



dt 



£dr=- 



di 



c (k ■ p) 



dRQED = huj'dWfi, 



{oj'fdu'dn _ _ uj'd 2 k' ± d(k ■ k') 



c(k ■ k') 

one can reduce Eq. (1321) for the polarized part of the emis- 
sion to the same form as that of Eq. (Q2 



dui'dn 4ir 2 c 



(Aqed(w') • I*) | , 



where 



Aqed(^') 



* + v(t) I w . 
exp < — — [t 



C 



fi 



(n.r(t)). 



C fi 



dt. 



where t— , t + are the time instants when the electron 
crosses the hypersurfaces £ = £_ and £ = £ + correspond- 
ingly. In the considered strong field case the finite inte- 
gration limits are admissible as long as the integration 
span well exceeds the formation time. Therefore, 



fi 




Aqed(w') = W-^-Aci — 
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^qed(^) n dUu'/C fl ) 



duj'dn 



fi- 



duj'dn 



(33) 



Thus, the QED effect on the emission from an electron 
reduces to the classical electric-dipole emission from a 
moving charge in an electromagnetic field with a very 
simple rule to transform the polarized emission intensity 
and polarized emission amplitude, which accounts for the 
recoil effect. 

However, the electron in the QED-strong field emits 
not only as an electric charge, because it also possesses a 
magnetic moment associated with its spin. In Eq. (|32l) a 
depolarized contribution to the emission is present. This 
contribution may be related to the electron spin, which 
is assumed to be depolarized. The depolarized emission 
energy related to the interval of frequency and to the 
element of solid angle and summed over two polarization 
directions (i.e., multiplied by two) equals: 



.ndcpol 
QED 

du'dn 



(q/) 2 (l-C /< r 

47T 2 C 



2C* 



fi 



]T|(A QED ( w ')-r)| 2 + ^QED| 



<y> —4 




-2 -1 

log, (>/ " c ) 



FIG. 2: Emission spectra for various values of \- 




-2 2 

i°g 10(4/4) 



whcr 



[IF f*+ f ioj' , (n.r(()),l dt 

is a scalar amplitude of emission, introduced in a way 
similar to that for the vector amplitude. After derivations 
analogous to those of Section II, the radiation loss rate 
due to the electron magnetic moment reads: 



dl 



dcpol 
QED 



4i(r) 9/3 



dui' 



8- 



(1 - C fi ) 



2 _ro_ 

Cfi 



K 



2/3 



Cfi 



(34) 



Thus, the QED effect in the emission from an electron 
in a strong electromagnetic field reduces to a downshift 
in frequency accompanied by an extra contribution from 
the magnetic moment of electron. Note a general char- 
acter of these conclusions: only in the recoil parameter, 
Cfi, is there a direct dependence on the ID wave vec- 
tor. This dependence can be also excluded, because, for 
the photons emitted along the direction of the particle 
motion, the following approximation is valid: 



hoj' 

£m e c 2 



1 - Xro, 



Cfi 



'0 



1 



X r o 



so that in QED-strong fields the emission spectrum is 
also universal and may be parameterized with the sole 
parameter, I c \. 

Combining Eqs. ()19I30I22[) one can derive the condition 



- £- » 



|da/del' 



(35) 



FIG. 3: Emitted radiation power in the QED approach vs 
classical (solid); an interpolation formula Jqed = J c i/(1 + 
1.04 v / 7 c i/J c ) 4/3 (dashed). 



Under this condition, the time interval within which the 
emitting electron locates between the wave fronts £ = £_ 
and £ = £+ much exceeds the formation time of emission 
validating the above considerations. 

With these results the scheme to account for the high- 
frequency emission as outlined in Section II may be eas- 
ily extended to 3D QED-strong fields. However, the ra- 
diation back reaction needs to be incorporated for this 
scheme to be consistent. 



IV. RADIATION AND ITS BACK-REACTION 

Unless the field is QED-strong the radiation back- 
reaction in a relativistically strong laser wave may or 
may not be significant. The condition for the field to 
be radiation-dominant (see, e.g., Ref.[H) is formulated in 
terms of the ratio between the magnitudes of the Lorentz 
force and the radiation force, which become comparable 
at intensities J ~ 10 23 W/cm 2 . For an electron mov- 
ing toward the laser wave, the radiation force starts to 
dominate at a lower wave intensity, depending on the 
electron energy H . The radiation back- reaction deceler- 
ates such an electron, the effect being more pronounced 
for longer laser pulses [14] . As the result, at intensi- 
ties J ~ 10 22 W/cm 2 the radiation back-reaction drasti- 
cally changes the character of the laser pulse interaction 
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with dense plasmas and 7-ray emission becomes a leading 
mechanism of the laser energy losses [25( . 

In QED-strong fields the radiation back-reaction is al- 
ways significant, as long as at each photon emission the 
electron looses a noticeable fraction of its momentum and 
energy. The matter of principle is also a consistency of 
the perturbation theory of emission in QED-strong fields. 
Within the framework of classical theory the momentum- 
energy change resulting from the radiation back-reaction 
should be small in some sense, to properly approximate 
the radiation force (see [3,[Ii| as well as the considera- 
tions relating to the estimate as in Eq. (PmO . In QED- 
strong fields this change cannot be claimed to be small, 
but the probability of emission can be! Specifically, the 
difference, — should be small enough, so that the 
probability of emission within the time interval of Eq. (|30[) 
should be much less (or at least less) than unity: 



Under these bounds, first, the condition (|35|) is satisfied. 
Therefore, the time interval ()30j) is much greater than 
the formation time and the emission probability is linear 
in 

dw fi (Z-,z + ) = (dw/dO(z + -t-). 

Second, the emission probability satisfies the condition 
(f3"6"| . Therefore, perturbation theory is applicable. In 
addition, the emission probability can be expressed in 
terms of the local electric field. Note, that consistency 
in (|37j) is ensured in relativistically strong electromag- 
netic fields as long as a < 1, with no restriction on the 
magnitude of the electromagnetic field experienced by an 
electron. 

Under the condition (|37p the probability may be ex- 
pressed in terms of McDonald functions: 

dW fl _ a * (/r~ K i (y) d y + r o r xX 2 Ki (r x )) 



A. Emission probability and radiation loss rate 

The derivations performed in Section II for the radia- 
tion loss rate, namely, the approximation of the angular 
distribution with the Dirac function and approximation 
of the frequency spectrum with the MacDonald functions 
may be directly applied to the emission probability. 

On developing the dot-product, {k' -pi), in T(£) in the 
TST metric, G^ v % one can find: 



T(0T(-6)=exp [i(Ti+T 2 )] 



where 



T, = 



(k-Pi) 



(k ■ k') 



2X c (k ■ k')(k ■ pf) \(k- Pi ) 

(fc-fc'){(e-a) + ^[a(6)-(a)] 2 rf6} 



T 2 = 



2X c (k-pi)(k-p f ) 



a = 



P±i) = P-LOi - (a) 



e-6 ' 

Integration over d 2 ]s.' ± then gives: 

dW fi {£-,Z + ) -It Jt^^^^i)^! 
d{k-k') 2wXc(k ■ Pi ) 2 

In strong fields the following estimates may be applied: 



(k ■ k') ~ X c {k ■ pif 



da 

di 



dWfi 



da 



Now the bounds for £ + — £_ can be consistently intro- 
duced: 

Ida/t^r 1 < £+ - £_ < min (V 1 (da/^l" 1 , l) • (37) 



(k ■ k') 
X(k-Pi) 



\i k -Pi) 



da 



An = \ ~z 



Ic 



Probability (similar to that found in [la]), is expressed 
in terms of functions of ro, and related to interval of dro. 
The way to introduce ro and x looks different from that 
adopted in Eqs. (117ll8[) . however, the difference is negli- 
gible as long as uj'/Si ~ (k ■ k')/{k ■ pi) for collinear k' 
and pi. The momentum of the emitted radiation, related 
to the interval of the electron proper time, may be found 
from Eqs.j30l38l): 



d Plad = f k , <k-m)dw fi d(k k , )d2 ^ = 



dr 



d(k ■ k')d 2 k ± d£, 



[p + k Odk-p,)- 1 )] / c(k'-k) 



dW 
' dr{)d^ 



drn 



(39) 



As with other 4-momenta, p ra d is related to m e c. To 
prove the 4- vector relationship (|3"9"|) . its components in 
the TST metric should be integrated over kj_ using the 
symmetry of T\. The small term, 0(l/(kpi)), arises 
from the electron rest mass energy and from the small 
(~ 1/|p_l|) but finite width of the photon angular distri- 
bution. In neglecting this term: 



dp, 



'rad 



dT 



dL 



QED 



drn 



dm 



'QED 



1 '<k-k')^dr 



being the radiation loss rate. The photon energy spec- 
trum, dlQw/dro, is described as a function only of the 
random scalar, ro, using only the parameter, x (see 
FigJ2]). The latter may be parameterized in terms of the 
radiation loss rate, evaluated within the framework of 
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FIG. 4: The emission spectrum for 600 MeV electrons inter- 
acting with 30- fs laser pulses of intensity 2 ■ W 22 W/cm 2 : with 
(solid) or without (dashed) accounting for the QED effects. 
Here hu c0 « 1.1 MeV for A = O.fyim. 



classical theory (see Eq.Q and FigEJ . The expressions 
for q{I c \) = Iqed/IcI and for the normalized spectrum 
function, Q(ro,x)> coincide with formulae known from 
the gyrosynchrotron emission theory (see §90 in [10]): 



9^3 
8tt 



dr r 



Ks (y)dy + r r xX 2 Kg (r x ) ], 



Q(ro, X) = 87rg ° \J (y)dy + r a r xX 2 K^ (r x ) 
As mentioned above, 

, j-pol , jdepol 

"-'QED _ a QED QED 



dr 



dm 



drn 



where polarized and non-polarized contributions are 
given by Eqs.([31[31 



B. Radiation back-reaction: radiation force 
approximation 

While emitting a photon, an electron also acquires 4- 
momentum from the external field (see Eq . (|3"TT) ) : 



(k ■ Pl ) - (k ■ k') 



k". 



The interaction with the field ensures that the total ef- 
fect of emission on the electron not to break the entity 
(j>f • Pf) = 1. As long as the angular distribution of 
emission is approximated with the Dirac function, the 
expression for dp F needs to be corrected to ensure ex- 
act momentum-energy conservation with approximated 
momentum of radiation. The choices of near-unity cor- 
rection coefficients in dpp are somewhat different in the 
cases x < 1 and \ 1. 



For moderate values of \ < 1 the radiation force, 
(dpF — dp ra .i)/dT, may be introduced. In this approx- 
imation it is admitted that the change in the electron 
momentum within the infinitesimal time interval is also 
infinitesimal. This 'Newton's law' approximation is per- 
tinent to classical physics and it both ignores the point 
that the change in the electron momentum at x ^ 1 is 
essentially finite because of the finite momentum of emit- 
ted photon and breaks the low bound on the time interval 
presented in (I35[) . The approximation, however, is highly 
efficient and allows one to avoid time-consuming statis- 
tical simulations. The approximation error tends to zero 
as x — > 0, however, it is not huge at x ~ 1 and even at 
X = 10. The latter can be seen from Fig f5] given below 
in which the average relative change in the electron en- 
ergy in the course of single-photon emission is presented 
(assumed to be negligible within the radiation force ap- 
proximation). 

Within the radiation force approximation the best cor- 
rection is dp F w k^(k' ■ pi)/(k ■ Pi). The total radiation 
force may now be found by integrating dpp over d{k' ■ k): 



dr 



{k-Pi 



-'QED 

m P c 2 



(40) 



The radiation force maintains the abovementioned entity 
as long as (pi ■ d(pj — p.^/dr) — 0. Eq. ptj)) is presented 
in a form which is applicable both with dimensionless or 
with dimensional momenta. 

In [13, [3 it was mentioned, that QED is not compat- 
ible with the traditional approach to the radiation force 
in classical electrodynamics and an alternative equation 
of motion for a radiating electron was suggested: 



& -or* /QED 



P^+ro^wn^pP, (4i) 



dr m e c 

where = eF^" /(m e c), F^ v is the field tensor and 

t q = 2e 2 /(3m e c 3 ) = (2/3)aA c /c. 

In the ID plane wave the particular expression for the 
radiation force can be found using the following equation: 

r ^^/ = k» (P 2 ^ Jc1 , . 

m e c z {k ■ p) 

With this account, the radiation force in Eq. f4"T|) is the 
same as its QED formulation in Eq. (|40| . This proves that 
the earlier derived Eg. ATI) has a wide range of applica- 
bility including electron quasi-classical motion in QED- 
strong fields and in the particular case of ID wave fields 
it can be directly derived from the QED principles. Note, 
that the efforts to derive the radiation force from quan- 
tum mechanics were applied many times (see [261 ] . the 
most convincing approach which give s the equation quite 
similar to (l4"TT) may be found in [27j|). However, for the 
first time the derivation from QED side is provided, with 
the resulting equation being different from those given in 
textbooks 
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FIG. 5: Expectation of the emitted photon energy. In dimen- 
sional units AS = hu}' and £ is the dimensional energy of an 
electron prior to emission. 



The way to solve Eq.(|4Tj) within the PIC scheme and 
integrate the emission is described in [l4j]. In Fig 0] we 
show the numerical result for an electron interacting with 
a laser pulse. We see that the QED effects essentially 
modify the radiation spectrum even with laser intensities 
which are already achieved. 



C. Radiation back-reaction: Monte-Carlo approach 

The radiation force approximation does not fully ac- 
count for the statistical character of the emission pro- 
cess at % > 1. Specifically, we mentioned above that 
in the 'Newton's law' approximation, the force, f , pro- 
vides only the infinitesimal change in the electron mo- 
mentum, Ap = fAi —> over an infinitely short time 
interval, as At —> 0. For radiation processes in QED- 
strong fields this point is in contradiction with a small 
probability, At ■ dWfi/dt — > 0, for an electron to acquire 
a finite change in momentum, \5p\ ~ |p|, in the course of 
emission. 

A more quantitative, though more cumbersome, de- 
scription may be achieved within the QED Monte-Carlo 
approach. It is convenient to relate the emission proba- 
bility to an interval of proper time, At = At /Si. From 
it follows that: 



dWfj _ /qed Q(rp,x) 
dr Q dr m e c 2 x r o 



(42) 



(Note, that on multiplying Eq. (|42l) by m e c 2 (uj' /£i) sa 
m e c 2 x r Qi one obtains again the formula for spectral dis- 
tribution of energy emitted per interval of time.) As long 
as Eq. (|4"2"f for the differential probability is available, one 
can find the expected energy of the emitted photon: 

J_ ,n = X 

£i ^ Jo 1/X ^o,x)(*oA-o)' 

A plot of the expectation for the ratio, < oj' /Si >, vs x 
is given in FigjSl with energy of emitted photons being 



denoted as AS. The total probability of emission within 
the interval of proper time is given by a complete integral 
of probability: 

W fi = Ar r^drdr = Ar I -^( ( fY\ 
J ar m e & \t% / 

Both within the QED perturbation theory and within the 
Monte-Carlo scheme W is assumed to be small W < 1. 
The probability of no emission equals 1 — W > 0. The 
partial probability, Wfi{ui' < wj), for the emission with 
the photon energy not exceeding the given value, u' , is 
given by the incomplete probability integral: 



W fi {j 



"''^ dr 

Q{ro,X) • 

X r o 



Therefore, for a given interval of proper time and calcu- 
lated x, /S and Wfi < 1, the expression of the only 
scalar to gamble, co'/Si, in terms of a random number, 
< rnd < 1, is implicitly given by an integral equation 
as follows: 



, dr rnd Id 
Jo X r o Wfi \Si 



if the gambled value of rnd does not exceed W^: < 
rnd < Wfi. Otherwise, i.e. if Wfi < rnd < 1, the 
emission within this interval does not occur. 

Once the value of u>'/£ is found, the change in the elec- 
tron 4-momentum due to single photon emission during 
the time interval, At, may be determined as follows: 



a 
Pf 



(pi-Pi)[l-o//(2g)] 
{k- Pi )[l-uj'/S] 



P) 



It is easy to see that the identity (pf ■ Pf) = 1 is main- 
tained. To achieve this, a correction factor, 1 — u' / (2S) is 
applied to the momentum exchange with the wave field 
as present in Eq. (f3Tj) . 

The implementation of this method for 3D realistic 
laser fields together with simulation results and an ac- 
count for pair production will be described in detail in a 
forthcoming publication. 



V. CONCLUSION 

QED-strong fields in the focus of an ultra-bright laser 
may be realized, if desired, using the technologies which 
already exist. In any case, these effects will come into 
power when laser-plasma interactions are explored with 
the next generation of lasers. 

It is demonstrated that electron motion in very strong 
laser fields with pronounced QED effects may be success- 
fully described within the radiation force approximation. 
The necessary corrections in the radiation force and the 
emission spectra to account for the QED effects are pa- 
rameterized by the sole parameter, I c \. 
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